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FOREWORD 


Flame spread over solid combustible surfaces is a phenomenon of 
fiandamental importance with a significant application to fire safety. 

From a practical point of view, it is importemt to be able to predict 
the flame spread rate over a variety of surfaces of arbitrary config- 
uration and orientation and a range of environmental conditions. One 
of the aspects of flame propagation which is imperfectly understood is 
the role played by buoyancy in helping set the flame spread rate and 
its effect on extinction. There have been studies reported in which the 
buoyancy level has been varied by changing the ambient pressture and/or 
the inclination of the specimen. In one study, the results of experi- 
ments at' nearly zero gravity level have been reported. 

In work we report here, the buoyancy level has been changed by 
elevating the gravity level at which the experiments were conducted. The 
report describes the experimental setup and the results of a number of 
experiments conducted with it. For vertical downward flame spread over 
thermally thin paper samples, a dimensionless spread rate is found to 

It 

correlate with a Damkohler number which varies with the ratio of the 
pressure to the gravity level. Different correlations are found for dif- 
ferent ambient oxygen concentrations? for air the flames are extinguished 

when the Damkohler reaches a value of about 2 x 10 . This corresponds to 

2 

a pressure to gravity level of 2 kPa *s /m. 

A paper based on this work was presented at the Fall Meeting of the 
Western States Section/The Combustion Institute. It is referenced herein 

- i - 



as [21]. 


This report, except for a few pages replaced by this foreword, is 
identical to the senior authors' M.S.E. thesis in Mechanical Engineering. 
It represents an interim report on work sponsored under NASA-Lewis Grant 
NSG-3114 with T.H. Cochran as Project Monitor. Paul C. Shang is pleased 
to acknowledge partial support from an HEW Domestic Mining and Minerals 
Fellowship throughout the period of his studies and research. D.C. Win- 
chester and P. Fluty were helpful collaborators in collecting the exper- 
imental data. Ms. Carol Marsh typed several versions of the report in- 
cluding the final copy. 
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CHAPTER I 
INTRODUCTION 

A, Brief Background on the Flame Spreading Problem 

Flame spread over the surface of solid combustibles is of interest 
with respect to fire safety. Thus, a substantial amoxant of cinalytical 
and experimental work has been done in the past decade [1,2] *. However^ 
due to the complexity of the flame spread phenomenon, complete under- 
standing has hot yet been achieved. The purpose of this investigation, 
which is mainly experimental in nature, is to provide information which 
will aid in the Tinder standing of how buoyancy affects the flame spread 
phenomenon. 

Figure 1 shows a schematic diagram of the problem associated with 
this work. A fuel bed is ignited at one edge, and following an initial 
transient period, a flame propagates steadily down the bed in a two 
dimensional fashion. Heat released by the flame and transferred for- 
ward of it vaporizes virgin fuel, producing the fuel vapor necessary to 
sustain the flame. In Fig. 1, a flow, which may be either naturally in- 
duced or forced, opposes the spreading flame. 

The mechanism by which the flame in Fig. 1 propagates is by trans- 
ferring heat to the unburnt fuel. For thin cellulosic specimens burnt 

* Numbers in square brackets refer to literature citations in the Refer 
ence Section. Reviews by Sirignano [1] and Williams [2] give excel- 
lent discussions on the majority of those works. Only material 
pertinent to the present study will be discussed here. 
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vertically downward, heat transfer occurs by conduction in the gas 
phase [3] . Since no forward convective heat transfer is present, the 
opposing flow must be such that a local Peclet number is of order unity. 

A fuel bed may not be heated appreciably prior to flame arrival. 

In this case, the bed is thermally thick. Lastrina et al. [4] charac- 
terized the fuel bed thickness by a dimensionless group, 

T = (1.1) 

where T is the dimensional bed half thickness, 6 is the length ahead of 
the flame in which the temperature of the fuel bed surface rises from 
its initial temperature to a "burning” temperature, a is the thermal 
diffusivity of the fuel bed and is the flame spread rate. A fuel 
bed is classified as thermally thick if T > 1. Otherwise, the bed is 
thermally thin. Since the ignition length 6 is dependent on the gas 
phase oxygen mole fraction and pressure, a particular fuel may be either 
theamially thick or thin depending on the environmental parameters. As 
a result, the usefulness of (1.1) is limited by the presence of this 
artificial ignition length. 

Since pressure, fuel bed thickness, oxygen mole fraction and an op- 
posing flow influence the flame spread rate, a great deal of attention 
has been given to the problem of determining the effects of these 
variables on the spread rate. 



-4- 


McAlevy aind Magee [5] and Lastrina et al. [4] have found that flame 
spread rate data can be correlated with the empirical equation 

( 1 . 2 ) 

where P is the environment pressure, is the environment oxygen mass 
fraction, and m and ^ are both dependent upon the type of specimen and 
diluent gas. For thermally thin cellulosic specimens burnt vertically 
downward, they reported a weak pressure dependence ($-0.1) and a near 
unity mass fraction dependence (m$ 1) . 

Frey and T'ien [6] studied the effects of pressure, sample size and 
oxygen mole fraction using thermally thin paper samples. They observed 
that power law correlations such as (1.2) were inadequate for flames 
spreading near extinction. 

The effects of an opposed flow on the flame spread rate for thin 
cellulosic specimens have been studied experimentally by Sibulkin et 
al. [7] and Hirano and Sato [8] . Sibulkin et al. obtained their data in 
a horizontal wind tunnel. Their results show that spread rate is nearly 
independent of free stream velocity, U^, for 15 cm/s < < 55 cm/s. 

As the free stream velocity is increased further, gradually decreases 
until, for ^ 125 cm/s, normal two dimensional burning no longer occurs. 
Data obtained by Hirano and Sato in a vertical wind txinnel show the same 
type of behavior. 

The effects of buoyancy on flame spread over thin cellulosic speci- 
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mens have been studied by inclining the . burning sample with respect to 
the vertical [7], thus generating a component of earth normal gravity 
along the direction of spread whose magnitude varies with inclination 
angle. However, one must interpret these data with caution as the usual 
fluid mechanical symmetry with respect to the specimen axis of symmetry 
is no longer present. It has been observed that at large angles to the 
vertical (usually > 30®) the flame on the underside of the specimen be- 
comes unstable, and fluid flow there may be in the same direction as that 
of flame propagation 19, 10] . 

McAlevy and Magee [5] have shown that for thin cellulosic specimens, 
the vertical downward flame spread rate is inversely proportional to the 
specimen thickness. They obtained specimens of varying thickness by 
fusing wetted index cards together under high pressure. However, Fernan- 
dez-Pello and Williams [ll] reported that two sheets of cellulosic material 
bonded together would spread apart under heating so that each sheet burnt 
independent of the other. They also reported that laminated sheets formed 
from, more than two individual sheets or a single sheet placed against an 
insulating backing did not exhibit self sustained downward propagation. 

Workers at NASA obtained spread rate data for thin cellulose acetate 
specimens burnt in zero and earth normal gravity environments [12, 13, 14]. 
In general, their results show that removal of buoyancy causes the spread 
rate to drop. However, since the usual product removal mechanism is ab- 
sent in zero gravity flames, one must be careful in comparing zero and 
earth normal gravity data. The ' buildup of combustion products may par- 
tially quench the flame, causing it to propagate more slowly [15] . 
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Scme analytical work pertinent to the btiming of thin cellulosic 
specimens is available. De Ris 116] assumed that the flame propagating 
over a fuel bed could be viewed as a thin diffusion flame. A uniform op- 
posing flow was considered. Consequently, no forward convective heat 
tremsfer was allowed. By neglectijig radiation effects and heat transfer 
in the fuel bed, the only forward heat transfer mechanism was gas phase 
conduction. A single global reaction was assumed to take place in the 
gas phase 



V* (Fuel) + V (Oxygen) (Product 1) 

f o . p 


(1.3) 


+ V" (Product 2) + HEAT 
P 


where v is the stoichiometric coefficient and terms enclosed by paren- 
theses are chemical formulas. Defining Schvab-Zeldovich vauricJsles, De 
Ris found the following expression for the flame spread rate; 


V xp C (T - T )= k(T^ - T ) 
F s ps vap f vap 


(1.4) 


where 


T = T + (BL/C ) (1 - lA) - (T - T ) 
f vap ~ P vap “ 


(1.5) 


K = [B/ln(l + B)]'^ Q/(MpV^L) 


( 1 . 6 ) 


Here is the fuel bed density, is the specific heat of the fuel bed, 

k is the gas phase thermal conductivity, is the fuel bed vaporiza- 

tion temperature, is the environment temperature , is the downstream 

asymptotic flame temperature, 0 is the heat of combustion of MV* grams of 

^ F f 

fuel, L is the latent heat of vaporization and _B is the mass transfer 



driving force. Note that in this expression V is independent of the op- 
posing flow velocity. More recently. Prey and T’ien H7, 18] considered 
the same problem. In their formulation, the gas phase momentum equa- 
tion was dispensed with through the use of the Oseen approximation. A 
one step second order reaction was assumed for gas phase oxidation, and 
a first order pyrolysis reaction was used for solid phase decomposition. 
Their numerical results show to decrease continuously with an increas- 
ing uniform opposing flow, V . Although the Solid and gas phase reaction 
rates and mechanisms are not Icnown exactly, Frley and T'ien were cUsle to 
predict flame extinction at low Damkohler numbers . The different depend- 
ence of V on the opposing flow velocity is because different chemical 
kinetics were used. De Ris assumed infinitely fast chemical kinetics; 

Prey and T'ien used finite rate chemical kinetics. 

B. Objectives of Present Study 

There is a sxibst«uitial amount of literature on flame spread over 
thin solid conibustibles . Most studies have shown experimental results and 
some correlations of the results. However, the effect of buoyancy has 
mostly been studied by inclining the fuel bed or by varying environmental 
pressTire in the normal gravity environment of the earth. The only study 
in which the effects of buoyancy were investigated by varying gravity 
level has been the NASA zero gravity work. 

The objectives of this investigation are to determine the effects of 
buoyancy on the flame spread rate over thin paper samples. Changes in 
pressure and gravity level affect the buoyant flow which is induced by 
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the hot downstream flame. The effects of these two parameters on the 
vertical downward flame spread rate over paper samples burnt in air and 
in a mixture which was 50% ©2/50% N 2 by volinne were studied. 

Previous correlations of flame spread resTilts have been limited 
to some ad hoc empirical equations. It is also an objective of this 
study to correlate flame spread rates and independent physicetl parameters 
in dimensionless groups. 
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axial cable for tran.siBitt,.iiig talevision signals, 2) a closed circuit tele 
vision system and a tape recorder, and 3} a reaote ignition and time 
spread aeasurement systaa. 

3. Speciaien Holder and Spread Rate MeasureBient. . The specimen hold- 
er used in this study consisted of two alimimjia plates {0.079 can thick). 
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was constantly purged by dry nitrogen from a compressed gas bottle. 

After the paper samples were dried, they were placed in an air tight 
plastic jar containing a molecular sieve and telltale silica gel as a 
desiccant so that a relative humidity of nearly 0% could be maintained. 
These techniques are widely used to provide uniform samples and are 
necessary. to achieve, reproducible results. 

The sample to be burnt was then placed between the two aluminum 
plates of the specimen holder. Three posts inside the chamber were used 
to mount, the., specimen holder , and its orientation was checked with a level 
prior to each run. 

A nichrome wire (gage 29, 13 fj/m) was used to ignite the paper 
sample. At elevated gravity, ignition of the paper sample using the 
igniter wire alone failed. Possibly, pyrolysis products were convected 
away by the buoyant flow induced by the igniter before ignition took 
place. In these instances a paper match was placed between the igniter 
and paper sample. The wire would first ignite the match head which would 
in turn ignite the sample. ' 

To fill the test chamber to the desired pressure with a preselected 
gas mixture, the chamber was sealed, and vacuum and gas supply systems 
were connected. To insure that contamination was less than 1%, the test 
chamber was evacuated and purged before it was filled to the desired 
pressure. Chamber pressure was monitored by appropriate test gauges 
manufactiired by Ametek. These pressure gauges were cleaned for oxygen 
service, and their accuracy was certified by the manufacturer to be 
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within 0.25% of full scale. 

For earth normal gravity experiments, the spread rates were obtained 
by the technique described above. For elevated gravity experiments, 
additional preparation was required. The supporting frame was locked 
at an appropriate angle of tilt. This angle, which is a function of 
gravity level, was calculated so that the resultant. of the earth normal 
gravity and the centrifugal acceleration was parallel to the s\irface 
of the paper sample for vertical downward flame spread. The actual tilt 
angle did not deviate from that calculated by more than 0.5 degrees. 

Gravity levels were controlled by varying the rotational speed of 

the centrifuge at a fixed radius.* A linear accelerometer (Schaevitz 

Model* LSB) whose output was amplified to give a 1 v/g reading on a 

e 

digital voltmeter was mounted on the frame to indicate gravity level. 

The accelerometer was calibrated directly by rotating it, thus produc- 
ing components of different magnitudes along its sensing axis. There 
is a maximum error of 0.3% in the accelerometer output. The accelero- 
meter has a ^ 25 g^ linear operating range. When the centrifuge had 
reached the desired rotational speed, the mechanical method described 
earlier was used to obtain the spread rate. 

The air used in the experiments was supplied by a compressor. 
Moisture in the compressed air was removed by a heatless dryer. The 50% 
0^/50% gas was supplied from a compressed gas cylinder whose contents 
were certified by Matheson Gas Products to be within 0.02% of the desired 
volume concentration. 



CHAPTER III 

RESULTS AND DISCUSSION 
A^ Results 

During the earth normal gravity experiments, several motion pic- 
tures were taken. They cover a pressure ranger of" 68.93 kPa; up to 413.58 
kPa for both index cards and paper tape burnt vertically downward in air. 
Motion pictures were also taken at varying gravity levels. Results of 
these motion pictures are discussed in the section on photographic ob- 
servations . 

About 90 vertical downward spread data for different gravities 
were recorded - - 58 in air, and 32 in the 50% 0^/50% N^ mixtxire. These 
data are tabulated both in raw and dimensionless forms in Appendix C. 

B. Power Law Correlation of Data 


The pressure dependence of normal gravity spread rate data is usual- 
ly correlated by power law relations [4, 5, 6] . Such results for verti- 
cal downward and horizontal spread from the present study are displayed in 
Figs. 8 through 11 to demonstrate that our measurement techniques are 
capable of reproducing results found in the literature. The straight 
lines shown in these figures are a calculated least square error fit to 
the data. 

The pressure exponent of about 0.1 for vertical downward spread in 
air (Figs. 8 and 9) is typical of values found in the literature [5,6]. 
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Figure 10. Pressure dependence of for index cards (T - 0.0098 cm) 
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The pressure dependence Is smaller for the thinner specimens. The 
smaller pressure escponent for data obtained in a 50% 0^/50% mix- 
ture is parobeQsly because conditions are further from extinction as sug- 
gested by Frey and T'ien [6] , 

The pressure dependence for horizontal spread (Figs. 10 cuid 11) is 
dependent on the specimen thickness as well as the oxygen mole fraction. 
In general, the pressure es^onent decreases with increasing oxygen mole 
fraction and decreasing specimen thickness for the same pressure range. 
Similar observations can be foxind in Ref. [6] . 

As noted previously, in the case of a paper sample,^ some investi- 
gators [9, 10] have observed unsteady fluid flow beneath the svurface in 
the same direction as flame propagation, causing convective heat trans- 
fer ahead of the propagating flame. This results in a much faster spread 
rate and larger pyrolysis zone than for vertical downward flame spread. 
Hirano et al. [9] observed that this instability is dependent on the con- 
figuration euid position of the remaining ash. Kashiwagi and Newman [10] 
reported that, for cx - cellulose sheets (38 cm x 84 cm) shallow cuts nor- 
mal to the direction of spread along the surface increased the frequency 
at vdiich the ash shcirds fell. The unstable lower flame motion was then 
not obseanred. In the horizontal flame spread experiments done in oiir 
study, curling of the ash did not affect the flame, and the data were re- 
producible. The gap width allowed for flame propagation was 1.92 cm 
for all data reported. When the gap width was increased to 10 cm, a 
width compcorable to that used by Hirano et al., curling of the ash stcirt- 
ed to affect the spread rate . For these wide samples , the remaining ash 
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divided into several portions with each portion assuming a different posi 
tion. When a portion of the ash hung below the specimen stirface^ more 
fluid was forced toward the xinburned material than when the ash was 
above the specimen stirface. Also, the configuration of the remaining 
ash was time dependent. These two effects together give rise to er- 
ratic burning and an uneven pyrolysis front. 

C. Dimensionless Correlation of Data 

Correlations such as those displayed in Figs. 3 through 11 were ob- 
tained by varying one parameter while other parameters that may affect 
the spread rate were held constant. Such correlations are not partic- 
ularly useful in attempting to determine the mechanism responsible for 
the behavior of the spreading flames. In this section, a method is pre- 
sented of correlating in dimensionless form the physical parameters that 
affect flame spread. 

In Figs . 12 through 14 , dimensionless correlations of the vertical 
downward flame spread rate are displayed. The dimensionless groups used 
here to correlate spread rates appear as dimensionless parameters in the 
mathematical formulation of the vertical downwaird flame spread problem 
[21] . In each figure the cibscissa is a Damkohler number; the ordinate is 
analogous to the fuel/air ratio in the region of flame attachment. The 
mass fraction of oxidizer also appecirs as a dimensionless parameter in 
the problem. Therefore, a different curve is generated for each mass 
fraction of oxygen. Values used for the various physical parameters are 
listed in Appendix C . 
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Figure 12. Dimensionless correlation of flame spread 
rate for index cards burnt vertically 
downward in air 
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Fi^re 13 . Dimensionless correlation of flame spread rate 
for adding machine tape burnt vertically 
downward in air 
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Figure 14. Dimensionless correlation of flame spread rate 
for index cards and adding machine tape burnt 
vertically downward in a 50% 0^/50% mixture 
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Frey cind T'ien 118, 19] have used a similar correlation. In their 

If 

theoretical ainalysis, they obtained a Damkohler number 


B = 


Bk 


CM V 
p ox a 


C3.1) 


as a dimensionless parameter in the energy equation. Here k and C are 

P 

the gas phase thermal conductivity and specific heat at constant pres- 
sure , B is the gas phase pre-exponential factor of a presumed second 
order reaction and is the molecular weight of oxygen. A uniform 
flow velocity, V , opposing the flame, was implied through the use of an 

©seen approximation. Their correlation was essentially a plot of V vs. 

F 

B for a particular fuel thickness and density. The correlation presented 
here is dimensionless and includes these variations. 


For the vertical downward flame spread case studied here, the fluid 
flow opposing the flame is a buoyant flow induced by the hot flame 
downstream of the point of flame attachment. The flow ahead of the flame 
must be sufficiently slow so that forward heat transfer by conduction 
through the. gas phase is capable of sustaining propagation. Consequently, 
the usual boundary layer approximation for natioral convection on a verti- 
cal flat plate is inapplicable. In the correlation presented here, the 

opposing flow velocity, V , is approximated by the characteristic buoyant 

a 

y. 

velocity, (VgBAT) . Here V is the gas phase kinematic viscosity, g is 
the acceleration of gravity, 3 is the volumetric coefficient of expansion 
and At is a characteristic temperature difference. This characteristic 
velocity is appropriate because in the vertical downward flame spread 
problem, as in the "general", vertical flat plate natural convection prob- 



-32- 


le?n discussed by Hellunjs. ^d Churchill 122 J , there is no . real characr* 
teristic length. The DarokShler huiqber defined in Eqn, 3,1 can therefore 
be written as , 


B = 


' • Bk • 

C M (VgBAT)% 
p ox 


(3.2) 


The Damkohler mimber is the ratio of a fluid mechanical time, t^, and 

the chemical time, t • For a unity Prandtl n\mber, the Damkohler mmiber 

c 

becomes 


1 / 

(V^/ggAT) ^ 


n = _ 

M (VggAT) yi 

ox 


(3.3) 


-1 y 

Here (Bp/M ) is the chemical time. The ratio of (V^/ggAT) , the 

OX 

buoyant length, and (VgBAT) , the buoyant velocity, represents the 
fluid mechanical time. The gas phase density is denoted by p. 


II 

As gravity increases or pressure decreases, the Damkohler number 

decreases • Since a decrease in Damkohler n\imber means an increase in 

the chemical with respect to the fluid mechanical time, reactants have 

time to diffuse through the reaction zone, causing the flame temperature 

to drop. As a result, less heat is available for conduction heat transfer 

to the virgin fuel. Therefore, if the initial transient is excluded, 

there exists a Damkohler number below which a steadily propagating two 

dimensional flame is not possible. For index cards and paper tape burnt 

vertically downward' in air, our results show that the flames are unable 

*» . 

to spread below a Damkohler number of about 2.0 x 10 . However, index 
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cards and paper tape extinguish at slightly different Dainkohler 
numbers. Possibly, this result can be an effect of slight differences 
in chemical composition which would cause different gas phase pre- 
exponential factors for the two cases. The same pre-exponential 
factor was used in correlating the data., Pigures 12 and 13 display 
such; results . As the. extinction limit was approached, one or more 
of the attempts to obtain a steadily propagating flame would fail. 

For those values of B, the results are shown as filled data points. 


The ordinate, p V t , is 

S F 


P V 

S F T 

P V T = — ^ "\y 

® p (VgBAT) (vVgBAT) 'V 


(3.4) 


or 


P 

— — — • S F 

DVT = 

^S F y 


(3.4a) 


where p x is the paper sample area density and y is the dynamic vis- 
s 


cosity. The niomerator, p V x/ is proportional to the amoiint of fuel 

s F 

available for combustion per unit sample width. The denominator, y , 
is proportional to the mass flux of gas mixture per unit sample width. 
Note that is independent of pressure and gravity level.. There- 

fore, the characteristic mass flux of gas mixture per unit width is 
independent of pressure and gravity level. It would be interesting to 
know whether or not the actual mass flux per unit sample width is also 
independent of pressure and gravity level. However, this would require 
complete knowledge of the gas phase velocity and temperature profiles. 
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As we showed in Fig. 6, for gap widths. ^ 1.92 cm, the flame spread 
rate for vertical downward flame spread in air up to three times earth 
normal gravity is xmaffected by specimen width when the sample is clamped 
in aluminxim plates. Also, Frey and T'ien [6] found the influence of gap 
width decreased as oxygen mole fraction was increased. The experimental 
results obtained here show that the numerical predictions of Frey and 
T'ien 118, 19] , which have V_ continuously decreasing as B decreases, 

r ■ 

are qualitatively correct. 

PgVpT as a function of B has a continuously decreasing slope with 

increasing B. As B , the paurameter p V T may become independent of 

s F 

B and possibly approach the thin flame limit of De Ris [16] . 


P„VT = ~ -. '^va p. (3.5) 

C T - T 
ps vap 

Using values of T and C from Frey and T'ien [18] and De Ris* down- 

stream asymptotic flame temperature, we found a value of 5.9 for air. 

Figures 12 and 13 show that p V T , as a function of B, appears to ap- 

s F 

proach an asymptotic value. If B is increased further, we expect the 
spread rate to become independent of the opposing buoyant flow. 

Portions of the present correlation for air that are away from the 

- ->9 -M 

extinction limit (B > 4 x lO** ) exhibit V <== g , or V “ V , Data from 

F Fa 

Sibulkin et al. [7] and Hiraino and Sato [8] show this type of dependence 

in certain ranges of V (Figs. 15 and 16) . In these wind tunnel experi- 

a 
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Figure 16, Dependence of Vp on the forced velocity V^. 

Data taken from Ref. [8], 
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ments, a laminar boiindary layer whose thickness varies as always 

present. From tiie data of Hirano and Sato [8] , the bomdary layer thick- 
ness appears to be substantially greater than the downstream flame stand- 
off distance. Consequently, for the lower velocity points, the flame 
would be expected to be entirely immersed in the slow flow region next 
to the surface. This observation may explain the independence of the 
flame spread rate on flow velocity for the lower free stream velocities 
in these experiments. The boundary layer thickness was calculated based 
on the sample length of Hirano and Sato. 


The V “ V relationship for the air results implies that, at a 
F a 

V 

constant gravity level, the flame spread rate is proportional to P'® . 


Such dependence is consistent with the pressure exponent of about 0.1 


reported by McAlevy and Magee [5] and with the data shown in Fig. 8. 

- 1/9 - 1 / 

Also, V oc g is implied by V a: V ' ^ . The inclined bed data of Sibul- 
F Fa 

kin et al. [7] exhibit this type of dependence at angles near vertical 


(Fig. 17). 



1 . Normal Gravity Photographs . Motion pictures of vertical down- 
ward flame spread on index cards and paper tape were taken at aeye.ral 
pressures. From these motion pictures, we saw that the flame color was 
greatly dependent on the pressure level. At high pressure (i.e. high 
Damkohler number) luminous flames occurred whose color was orange with a 
small bluish region near the area of flame attachment. The flame size 
was smaller than those at a lower pressure, indicating that the chemical 



0.2 0.3 0.4 0.5 0.6 C 

Sin e 


Figure 17. Dependence of on the orientation for thin 
paper samples. Data taken from Ref. [7]. 
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reactions were fast and only occurred in a small region. As the pressiare 
was reduced/ the residence time of the reactants in the flame zone was 
increased. The flame gradually transformed into a blue colored flame 
with a visible orange glowing char. This decrease in luminosity was 
caused by a local decrease in fuel/air ratio. Flame size was greatly 
enlarged at lower pressures. The flame shape in all cases was botindary 
layer like. 


2> Photographs with Identical Damkohler Number . Motion pictures 

of vertical downward flame spread on index cards were taken at B - 14 x 

10*^ and B-2.7 x 10** in an attempt to verify the present model . B - 14 x 

10** was obtained at 413.58 kPa and 4g / and at 101.33 kPa and Ig . 

^e e 

B - 2.7 X 10** was obtained at 68.93 kPa and 2g / and at 137.86 kPa and 

© 



These pictures revealed that flames with identical Damkohler 
mmiber did not have the same characteristics. Instead, flame size and 
color seemed to beiinaffected by changes in gravity level. This suggests 
that for flames of identical Damkohler number, only regions near the 
flame attachment area are identical. However, in our experiments , this 
region could not be photographed as the alimiinimi plates* thickness is 
about the same as the flame stand-off distance near the foot of flame 
attachment [9] . 



CHAPTER IV 


CONCLUSIONS 

The dependence of flame spread on buoyancy has been demonstrated 
by varying gravity level and pressure. The flame spread rate is shown 
to decrease as the buoyancy induced flow increases. 


A method of correlating flame spread data using dimensionless para- 
meters has been presented. Flame spread rate, V , is expressed by 


P 

S F 

p V T = — 


(4.1) 


which allows variations in sample thicJcness and density. The Damkohler 
momber 



Bk 

y: 

C M (Vg3AT) 
p ox 


(4.2) 


is the dependent variable. If the Damkohler number is reduced, the 
flame spread rate drops xintil extinction occiars"; The extinction limit 
is about B = 2.0 x 10*^ for paper samples burnt in air. 
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APPENDIX A 

ERROR IN EXPERIMENTS. 

!♦ Oxygen Mole Fraction Changes Druring a Run 

The most severe changes in oxygen mole fraction inside the chamber 
occur during the bviming of index cards in the 50% 02 / 50 % N 2 mixture 
at 20.68 kPa. For reasons explained in the Spread Rate Measurement 
section, a 7.52 cm by 1.92 cm strip of index card is considered here 
for analy^s ptirposes. 

Such a paper sample has a mass of 0.21 g (Appendix C) . Assuming 
pure cellulose, (^6®io*^5^x' following reaction to estimate 

the oxygen mole fraction changes: 

=6"l0°5 * ®°2 * ^“2 * ^»2° 

The moleculax weight of the cellulose monomer is 162 g/mol. Therefore, 

“3 — 3 

1,30 X 10 moles of cellulose and 7.82 x 10 moles of oxygen were con- 

s\imed during each experiment. At standard conditions, the test chamber 

contains 0.346 moles of gas. Hence, for a 50% 02/50% N 2 gas mixture at 

20.68 kPa, there are 0.173 moles of oxygen available for combustion. 

It follows that the total percentage change in the oxygen mole fraction 

is no more than 0.8%. 

2. Error Analysis 


The error in determining the flame spread rate, V , is the error 

•t 
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in spread rate measurement which includes two variables, t and d. The 
relation between these two variables can be written in the following 
form 



(A-1) 


where d and t are the distance and the time for the flame to travel the 

distance. Therefore, the error of V_ by measurement is the combination 

£ 

of the errors caused by measuring d and t. 

The variation of by measurement can be analyzed by the standard 

error analysis method [23] . Equation A-1 is differentiated and divided 

bv V„. Therefore, 

F 


AV -At Ad 


Taking the root mean square of the terms in Eqn, A-2, we, have 


AV, 


F^rms 



K 


(A- 3) 


where A V„ is the probaible error in V_. 

F,rms F 


The variation of t and d are discussed individually as follows: 

At 

— : The possible error in time measurement can arise from the following 

sources 

(1) The error in measuring the time. In hand timing the spread rate. 
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a reflex time of 0.1 s is assmed for each start of stop of the 
digital clock. The same amoiint of error is assmed for the re:- 
sponse of the microswitches used. The error is not more than 3%. 

(2) Uncertainty in reading clocks is less than 0.5%. 

(3) The error in the clocks is not more thcin 1%. 

^ : This error can only be caused by improper reading of the rules 

and by the error of the irulers themselves. As with the clocks, 
this is no more than [1^ + (0.5)^] =1.1%. 

From the above analysis, it is concluded that the probable error 
for determining V_, as expressed by equation A-3, is 

AV_ . 

- . = [3^ + (0.5)^ + + (l.l)^]/2% = 3.4% (A-4) 

F 

With the present correlation scheme, in which P«V X depends on B, 
errors of p V„x and B should also be discussed. In these two dimen- 
sionless parameters, errors in parameters that have been taken as constant 
are assumed to be zero. These parameters are the gas phase thermal con- 
ductivity, k, the gas phase dynamic viscosity , y , the gas phase specif- 
ics heat at constant pressure, C , the gas phase pre-exponential factor, 

XT 

B, the coefficient of expansion, 3 , and the characteristic temperature 
difference, At. 

The expression for the dimensionless group Pg^pT is repeated here 
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for convenience* 


P 

s F 


P == 

s F y 


(A-5) 


Similar to the error analysis for V^, the error of 


written as 


Ap^v^T Ap^t AV^ Ay 


(A-6) 


P V„T P T 

^s F s 


y 


Under careful calculation of properties used, it is reasonable to assimie 
that the error in the solid phase area density, p^T, is no more than 4%. 
Thus, the probable error of p V t, denoted by p V„x /p V_x,can be 
determined as follows 


P 

s F rms 


P 

S F 


r . \ 

Ap X 




2 - 

+ 


V 


v_ 


21 1 


% 








Ot 


= (4' 


+ 3.4^)^^% = 5. 


25% 


(A-7) 


The above value represents the maximimi probable error of the ordinate 
in the present correlation. 


The Damkohler number, B, as defined in Eqn. 3*2 is Bk/(C 

p ox 

zL — 

(vg$AT) / ) . As in the analysis above, the probable error of B can be 

determined by the following equation 
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(A-8) 


The environment pressure, P, is known to within 1%. The output 
of the digital voltmeter used to determine the gravity level under which 
the experiments were conducted showed an oscillation in gravity level 
due to an imperfection in regulating the rotational speed of the cen- 
trifuge. The deviation in gravity level, g, was found to be less than 
1.5%. Therefore, the maximum probable error of B, the abscissa of the 
dimensionless correlation developed here, is 


AB 

rms 

B- 




'2 

2 

' 2 



- • 1 

+ 

- • 1,5 



I3 J 


3 j 


^2 


% = 1 . 2 % 


(A- 9) 



APPENDIX B 


Effect of Coriolis Acceleration on Flame Spread Rates 

A Coriolis acceleration, 2o)x f, is always present for an object 
moving with a radial velocity f in a system rotating at an angular 
speed In all elevated gravity experiments described here, the test 
chamber was locked at an angle so that the resultant of the centrifugal 
acceleration and earth normal acceleration was in the same direction as 
the flame spread. Therefore, gas particles of the buoyant flow induced 
by the flame downstream of flame attachment have radial velocities, and 
thus e^erience a Coriolis acceleration. 

For the data reported here, the paper sample orientation was such 
that the Coriolis acceleration was normal to the sample surface. The 
magnitude of the Coriolis acceleration is estimated at 20) (Vg3AT) 

y 

which after introducing the relationship g = + g^] ' to eliminate 

0 ), becomes 2 [ (g^- g^)/R^] (Vg3AT) . At a radius of 6.1 m,the radius used 
to obtain the data shown in Figs. 12 through 14 , the ratio of this ac- 
celeration to the acceleration of gravity along the sample surface is 
0.06 for the worst case (i.e. at 4g and 20.68 kPa) . 

The effect of Coriolis acceleration on the results reported here 
was tested by varying the radius of rotation at constant gravity level. 
Figure Bl shows that, in general, the spread Vp is independent of Cori- 
olis acceleration , 

A Coriolis acceleration tangent to the sample surface can be ob- 
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Figure B1 Dependence of V on the Coriolis acceleration (index cards) 
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tained by rotating the chamber by 90 degrees about its axis of symmetry. 
Spread rate results (Appendix C) obtained with this configuration show 
a different extinction limit. The flame extinguishes at higher values 

It 

of the Damkohler number when the Coriolis acceleration vector is tan- 
gent, rather than normal, to the sample storface. Apparently, a cross 
flow,' due to the; tangential acceleration, blows the. flame out pre- 
maturely. 



APPENDIX C 

TABULATION AND REDUCTION OF DATA 
1. Data Used in correlations 

All gas phase physical properties were evaluated at 1500 K. 

BAT was set equal to 4. Values of the dynamic viscosity, y , the 
specific heat, C , and the thermal conductivity, k, of air were obtained 

Xr 

from Ref. [24] . Values of these physical properties for pure oxygen 
and pure nitrogen were obtained from Refs. [25] and [26]. The 50% © 2 / 
50% N 2 mixture's dynamic viscosity and the thermal conductivity were 
calculated by equations 9-5.4 and 10-6.9 of Ref. [27]. Molecular 
weights of air and 50% 02 / 50 % N 2 mixture are 28.97 g/mol and 30 g/mol, 
respectively. The ideal gas law was used throughout to calculate the 
gas phase density. 

The property values used are listed below: 

Air at 1500 K and atmospheric pressure: 

C ’= 1.230 x 10^ J/(kg*K) 

P , 

k = 0.0946 J/(s*m«K) 

-5 

]i = 5.40 X 10 kg/(m*s) 

50% 0^/50% mixture at 1500 K.and atmospheric pressure: 

Cp = 1.194 X 10^ J/(kg*K) 
k = 0.0931 J/(s*m»K) 

-5 

U = 5.63 X 10 kg/(m*s) 
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The area density,p^T/ of the paper tape and index cards was 
measured to be 2.890 x 10 ^g/qm^ and 7.381 x 10 ^gf/cm^, and the re- 
spective half thicknesses are 0.0043 cm and 0.0098 cm. The pre-ex- 
ponential factor, B = 1 X 10^2(303/ (mol * 3) , was taken from Ref. [19]. 

The vertical downward flame spread data have been arranged in order 
of increasing pressure for each gas mixture and s^ple thickness . Each 
data point listed represents an average of several spread measurements. 
An asterisk (*) indicates that one or more attempts yielded flames which 
could not propagate. 



-54- 


PressTire g/g V„ B 

(kPa) ® (cm/s) (x lo“) 


Air? Index Cards 


P T 
s F 


68.93 

1.00 

0.147 

4.301 

2.007 


1.50 

0.134 

3.282 

1.829 


1.75 

0.129 

2.962 

1.761 


2.00 

0.117 

2.709 

1.597 


2.50 

0.116* 

2.335 

1.584* 


3.00 

* 

2.068 

* 

86.16 

i.bo 

0.157 

4.991 

2,143 


1.50 

0.145 

3.809 

1.980 


2.00 

0.136 

3.144 

1.857 


2.50 

0.126* 

2.709 

■ 1,720* 


3. 00 

0.123* 

2.399 

1.679* 


3.50 

* 

2.165 

* 

101.33 

1.00 

0.154 

5.560 

2.102 


1.50 

0.145 

4.243 

1.980 


1.75 

0.142 

3.829 

1.939 


2.00 

0.134 

3.503 

1.829 


3.25 

0.122 

2.534 

1.666 


3.50 

0.115* 

2.412 

1.570* 


4.00 

0.112* 

2.207 

1.529* 


4.25 

* 

2.119 

* 

137'. 86 

1.00 

0.157 

6.827 

2.143 


2.00 

0.146 

4.301 

1.993 


3.00 

6,137 

3.282 

1.870 


4.00 

0.117 

2.709 

1.597 

275.72 

1.00 

0.177 

10.840 

2.416 


2.00 

0,157 

6.287 

2.143 


3.00 

0.145 

5.210 

1.980 


4,00 

0.140 

4.301 

1.911 
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413.58 

1.00 

2.00 

3.00 
3.50 

4.00 

0.178 

0.161 

0.154 

0.151 

0.145 

14,200 

8.946 

6.827 

6.160 

5.636 

2.430 

2.198 

2.102 

2.061 

1.980 

Air; Paper Tape 

68.93 

1.00 

0.369 

4.301 

1.985 


3.00 

0.283 

2.068 

1.522 ^ 


3.50 

0.280* ' 

1.866 

1.506* 


4.00 

* 

1.707 

* 

101.33 

1.00 

0.396 

5.560 

2.130 


1.50 

0.358 

4.243 

1.926 


2.00 

0.363 

3.503 

1.953 


2.50 

0.324 

3.019 

1.743 


3.00 

0.307. 

2.673 

1.651 


3.50 

0.298 

2.412 

1.603 


4.00 

0.285 

2.207 

1.533 

137.86 

1.00 

0.403 

6.827 

2.168 

275.72 

1.00 

0.419 

10.840 

2.254 

413.58 

• 1.00 

0.422 

14.200 

2.270 



50% 

02 / 50 % N 2 ; 

Index Cards 


28.68 

1.00 

0.525 

1.943 

6.995 


1.50 

0.499 

1.483 

6.649 


2.00 

0.469 

1.224 

6.249 


2.50 

0.460 

1.055 

6.129 


3.00 

0.447 

0.934 

5.956 


3.50 

0.424 

0.843 

5.650 


4.00 

0.468 

0.771 

6.236 

34.47 

1.00 

0.550 

2.732 

7.328 


1.50 

0.519 

2.085 

6.995 


2.00 

0.519 

1.721 

6.915 
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2.50 

0.512 

1.483 

6.822 


3.00 

0.512 

1.314 

6.822 


3.50 

0.492 

1.185 

6.549 


4.00 

0.493 

1.084 

6.569 

6 B .93 

1.00 

0.569 

4.337 

7.582 


2,00 

0.540 

2.732 

7.160 


3.00 

0.554 

2.085 

7.345 


4.00 

0.553 

1.721 

7 . 332 

101.33 

1.00 

0.602 

5.607 

7.982 


2.00 

0.600 

3.532 

7.955 


2.50 

0.566 

3.044 

7.542 


3.00 

0.547 

2.696 

7.284 


4.00 

0.539 

2.225 

7.177 

275.72 

1.00 

0.617 

10.930 

8.214 

551.58 

1.00 

0.646 

17.350 

8.601 


50% 

02 / 50 % N2 ,* Paper 

Tape 


20.68 

3.00 

1.221 

0.934 

6.495 

34.47 

1.00 

1.355 

2.732 

7.207 


3.00 

1.298 

1.313 

6.904 

68.93 

1.00 

1.430 

4.337 

7.606 


3.00 

1.324 

2.085 

7.042 

101.33 

1.00 

1.498 

5.607 

7.967 

275.72 

1.00 

1.502 

10.920 

7.987 

551.58 

1.00 

1.475 

17.350 

7.846 
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2, Elevated Gravity Data With 
Tangential Coriolis Acceleration 


The elevated gravity data for vertical downward flame spread with 
Coriolis acceleration tangent to the specimen surface are listed below 
in the order of increasing pressure. Values used for various parameters 
are the same as those used in reducing the data that appear in the pre- 
vious section. Again, these data are the average of several runs and 
asterisks indicate that one or more attempts to obtain steadily propa- 
gating flames failed. 


Pressure 

e 

(cm/s) 

B 

(x 10**) 

P= T 

"^S F 

Air ; Index , Cards 

68.93 

1.50 

* 

3.282 

* 

86.16 

1.50 

0.131 

3.809 

1.788 


1.75 

0.127 

3.437 

1.734 


2.00 

0.121 

3.144 

1.652 

101.33 

1.50 

0.145 

4.243 

1.980 


1.75 

0.141 

3.829 

1.925 


2.00 

0.129 

3.503 

1.761 


2,50 

0.126 

3.019 

1.720 


3.00 

0.113 

2.673 

1.543 


3.25 

* 

2.534 

* 







